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Abstract 

We prove a condition on / £ C 2 (K+,R) for the convexity of / o det on P§ym(n), namely 
that / o det is convex on PSym(n) if and only if 



*5" f" (s) + f'(s)>0 and /'(sXO VsGR+. 

' ' [ ns 



This generalizes the observation that C — In det C is convex as a function of C. 

1 Introduction 

The question of how to choose physically reasonable strain energy functions in nonlinear elasticity 
has attracted much attention and is not yet completely solved. The major breakthrough came 
with John Ball's seminal contributions [5J [TJ [3] introducing polyconvexity, i.e. convexity of the 
strain energy W as a function of the arguments (F, CofF, det F), see also [El [§]• Polyconvexity 
reconciles the physically reasonable growth condition W(F) — > oo as detF — > with the weak- 
lower-semicontinuity (quasiconvexity) , which in return implies ellipticity. A very simple example 
of a polyconvex function is the uni-constant compressible Neo-Hooke model 

W NH {F) = fJ,[(F T F -1,1) - 21ndet(F)], shear modulus [i > 0. 

The strain energy is isotropic, frame- indifferent, polyconvex, i.e. convex as a function of (i* 1 , det F), 
stress-free in the reference configuration and Wnh — > oo as det F — > 0. It is well known that the 
latter requirement excludes from the outset that F i-> Wmh{F) may be a convex function of F, 
[B]. Howewer, rewriting Wnh in terms of the Cauchy-Green deformation tensor C = F T F, which 
gives 

W NH (F) =W nh (C) =/x((C-l,l) - In det C) 

one may readily check that C i-> Wnh{C) is a convex function of C, despite its singularity in the 
determinant as det C — > 0. 

We surmise that convexity of the free energy w.r.t. C (or the stretch tensor U = V~C) is an 
additional, desirable feature of any free energy as it implies monotonicity of the stress-strain 
relation. 

In this short contribution we therefore investigate, which functions / g C 2 (IR+,]R) are such that 
C i— s- /(det C) is convex as function of C G P§ym(n) and generalize the well-known result that 
Cm- — In det C is convex on the set of positive definite symmetric matrices |101 HI [S] by proving: 
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Theorem 1.1. (A differential inequality characterization) Let f e C 2 (I 
the function 

/odet : PSym(n) — ► R , C i — ► /(detC) 

zs convex if and only if 



Then 



f"(s) + 



n — 1 



/'(s)> and /'(s)<0 Vse 



Proof. This is an immediate consequence of Lemmas 11.51 11.71 and 11.91 



(1) 



□ 



In the following we will reformulate the condition for convexity to obtain this result. We start 
with some preliminaries: 

By M nx ™ we denote the set of all n x n- matrices, §ym(n) stands for the set of all real symmetric 
n x n-matrices and PSym(n) for the set of all real symmetric positive definite n x n-matrices. 

Lemma 1.2. (Characterization of convexity) Let X be a normed space, g G C 2 (K, R) and 
K C X open and convex. Then 



g convex D 2 g(x).(z, z) Va; G K, z G span(i^) . 

Proof. [1, p.27] 

In particular we obtain 
Theorem 1.3. (Condition for convexity) For g G C 2 (PSym (n) , M) we have 
g convex <^=4> D 2 g(C).(H, H) ^ VC G PSym(n), iJ G §ym(n) 



(2) 
□ 



(3) 



Proof. Let -ftf := PSym(n) and X := §ym(n) in the previous lemma. PSym(n) is an open 
convex subset of the normed space §ym(n) (with operator norm): Use the characterization 
A G PSym(rj) (Ax , x) > OVx G R n \{0} and for convexity also the Cauchy-Schwarz- 

inequality. Furthermore span(i^) = span(PSym(n)) = §ym(n) = X: The inclusion „C " is obvious. 
For the other inclusion write A as a diagonal matrix (the corresponding transformation preserves 
positive definiteness and symmetry) and show that this can be written as a linear combination of 
positive definite symmetric matrices. □ 



By (A , B) = tr(AB T ) we denote the trace inner product of the matrices A and B. 



Theorem 1.4. (A condition for convexity) 
Let f G C 2 (R+,R). Then the function 



g := f o dct : PSym(n) 



C 



/(detC) 



is convex if and only if 



VC G PSym(n) WH G Sym(n) : 
/"(detC)detC+/'(detC)j (C" 1 , H) 2 - /'(det C)(HC- 1 , C~ l H) > 0. 



(4) 



Proof. Because / G C 2 ,det G C°°, also jeC 2 . It remains to be shown that 

D 2 g(C).(H, H) = det C ■ { [/"(dct C) • det C + /'(det C)] ■ (C^ 1 , H) 2 

-/'(det CO ■ (HC- 1 ,C" 1 H)} 



for C £ PSym(n) and H £ §ym(n), then the claim follows by Theorem ll.31 Because det is infinitely 
often differentiable on M nxn and (cf. [5]) Ddet{A).H = (Ad] A T ,H), where Adj A denotes the 
adjugate matrix of A. For invertible C and symmetric H we have D det(C).H — det C (C" 1 , ff ), 
and hence obtain by the chain rule 

Dg(C).H = Df {dot C)D dct(C).H = /'(det C) • det C • (C -1 , H) , 

and therefore, by chain rule and the fact that £>[C -1 ].iZ = —C HC , 

D 2 g(C).(H, H) = /"(det C) • (det C) 2 • (CT 1 , H) 2 + /'(det C) • det C • (C" 1 , H) 2 

+ /'(detC) -detC- (-C^HC' 1 ,H) (5) 
= /"(det C) • det C 2 ■ (C^ 1 , ff) 2 + /'(det C) • det C ■ (CT 1 , tf) 2 
- /'(det C) • det C • (iJCT 1 , C^ff ) . □ 

Lemma 1.5. The inequality f'(s) ^ Vs £ R + is necessary for Q. 

Proof. Assume s € R + satisfying f'(s) > 0. Let C — diag(l, 1, s) £ PSym(n) and H = 
diag(l, -1, 0, 0, ...) £ Sym(3). Then det C = s, (C" 1 , H) = 1 - 1 = and (HC- 1 , C~ l H) = 
(diag(l,-l,0,...,0) ,diag(l,-l,0, ...,0)) = 2. Together with (0} we obtain -2/'(s) ^ 0, a contra- 
diction to /'(s) > 0. □ 

Lemma 1.6. (j4]) /io/ds j/ and only if 



ViJeSym(n) VL^ 1 = diag(di, d n ) , 
where di, ...,<in G R+ <wid s -1 := det!? -1 = d\ ■ . . . ■ d n G R+ 

/"(«) + LW\ (D- 1 , H) 2 - ^ (D^H , if I?" 1 ) ^ . 



(6) 



Proof. Consider an arbitrary C £ PSym(n) in Then there is an orthogonal matrix Q, such 
that C = QDQ T C _1 = QD~ 1 Q T , where L> = diag(Ai, A„) and positive Aj. By the 

properties of the scalar product of matrices we have 

(C- 1 , H) = (QD~ l Q T , H) = (QD- 1 , HQ) = (D~ 1 ,Q T HQ). 

For H £ §ym(n) let H := Q T HQ and note that H varies over the whole of Sym(n) if and only if 
H does. Analogously, 

(HC- 1 , C~ l H) = (HQD~ l Q T , QD~ 1 Q T H) = (hD- 1 , D' 1 ^ = (d^H , ff-tr 1 } 
Denote di := A^ 1 and s := det C = det D = Yi7=i ^« an< ^ divide (0| by s > to obtain ©. □ 
Lemma 1.7. Let f £ C 2 (R + ,R) andfodet be convex on PSym(n). Thenf'(s) > Vs £ 



Proof. According to Lemma [1.61 © holds for all H £ Sym(n) and D 1 = diag(di, d n ). Let 
s € R+, fc e R\ {0} and iJ = k ■ D' 1 , as well as D^ 1 = diag(s - ~, 



<fc 2 ((/"( S ) + ^) (D- 1 , D- 1 ) 2 ^ (p- 1 ) 2 , p- 1 ) 2 )) 
=fc 2 ff/"( S ) + 



=fc 2 /"(S) + 



(ns- 2 /") 



s 

2 f 



/'(*) 



ns 4 ^'" I = n 



k 2 s- A ' n {nf'{s) + {n-l) f -^\ 



□ 
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For any matrix A let diag A be the matrix obtained from A by setting all non-diagonal entries 
zero. Let diag M „ X n be the set of all n x n— diagonal matrices. 

Lemma 1.8. For all P G diag M „ X n with non-negative entries only and all A € M nXTl the following 
holds: 

(P,A) = (P, diag A) =:a(P, A), (7) 
(PA,AP) ^ (P diag A, diag A P) =:a(P,A), (8) 
o- 2 (P,A) ^ n-a{P,A) . (9) 

Proof. Let P = diag(p x , ...,p n ), A = (ay)^- and calculate (P , A) = YTi=\Viaa = (P,diagA). 
Hence ([7]) holds. Direct calculation of PA and PA T yields 

n n 

(PA,AP)=tY(PAPA T )^Y,Py^+Y,T,P^ ka ^ > (P diag A, diag A P), 

i—1 i—1 k^i 

i.e. ©. For all P G diag M „ x „ and A G M" x ™, we have cr 2 (P, A) < n • ct(P,A). To see this, 
note that P and diag A commute, i.e. <r(P, A) = (P diag A , P diag A) = ||PdiagA|| 2 holds. By 
Cauchy-Schwarz-inequality this immediately implies 

cr 2 (P,A) = (P,diagA) 2 = (PdiagA,!) 2 ^ ||PdiagA|| 2 • ||1|| 2 =n-a(P,A) . □ 

Lemma 1.9. ((T|) is sufficient for the convexity of /odet. 

Proof. We will show To this end, let H G Sym(ri) and D^ 1 = diag(di, d„), where 

di, ...,d n G K+ and s := (d± ■ . .. ■ dn) -1 = detD arbitrary. Then P := D^ 1 and A :~ H satisfy 
all assumptions of the previous lemma. Using the notation from lemma [TTHl we can, without loss 
of generality, assume a(D^ 1 ,H) ^ 0, because otherwise (O becomes trivial by the assumption 
/' < 0. We denote a = o^D -1 , iT) and a = a{p~ x ,H) ^ (P _1 P , HD^ 1 ) by ©. Using ^ < 
and f"(s) + > by (TTJ) and 1 - < z=± we obtain @ from 

/'('A/n-1 ^ 2 /'OO /n-l, ^-IW^ , _2 /'(«) 



/" ( S ) + L^l (ij-i , _ H , PP- 1 ) > /» + 



CT • (7 



fj 2 



V n s ' 



2 Solutions to the differential inequalities 

In this section we are interested in the possible shape of the functions that satisfy ((1}. To make 
calculations and figures more concrete, we restrict ourselves to the case n = 3. 

Lemma 2.1. (Linear ODE) The linear initial value problem 

Ly := y' + g(x)y = 0, y(C) = n (LWP) 
on J — R + and where g{x) = ^ has one and only one solution. 

To find solutions to Lf > under the additional constraint y = /' < (which is equivalent 
to r] < because /' = is a solution) we consider the "limiting case": 

Lemma 2.2. (Limiting case for JTJ) TTie solutions to 

/!*(*) + I" • /L«(s) = and /L«(s)<0 V,sGR+ 
are given by /i imit : R+ — > R, s i— > c • s 1 / 3 + d, where c ^ , d G R. 
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Proof. Separation of variables gives the unique solution of (jLIVPp for £ > ^ 77: 



2 

3t 



dt J = 77 • cxp [ - - In - ] =77 £ 2/3 • x 



-2/3 



(10) 



Because 77 £ 2 / 3 sC 0, we have 7/ limit ^ 0, hence ([!]). The claim follows by integration of f' itait = J/Wm 
with c := 3?7£ 2//3 and constant d. □ 

If we consider an interval adjacent to £ on the left hand side, i.e. J := [£ — a, £], the conditions 
for a function y to be a sub- (or super) solution to y' — F(x,y), y(£) = 77 are 



F(x, v) in J , ^ 77 







^or It/ < 





F(x,w) in J, w(£) 77 respectively 



where (cf. (jLIVP]) ) F 2/3 (a;, y) := - ^ • y e C*°°(M + x K) yields 
Lemma 2.3. Le£ y 6e differentiable in R + , £ > ^ 77. XTien 



F 2 / 3 (x, y) , y(0 > 77 y(x) 



2/ llmit (x)=^ 2 / 3 .x- 2 / 3 



K,oo) 



Analogously: 



F 2 /3(x,y), ?/(£)^77 



f 2/limitW 0M ^ / 



By these considerations, we obtain information on the qualitative shape of the solutions to ([T]). 
(At first discussing the shape of y — /'.) Note that to fulfill y < 0, in Lemma |2~31 also ^ ^ 77 
must be satisfied. For 77 = 0, y limit = is the unique solution and intersects y in £. (0 > y(£) > 0) 
For y(£) > 77, y and y limit with initial value j/i im i t (0 = 2/(0 intersect in £. Hence we can 
consider t/(£) = 77 = j/ limit (£) only. Then 



implies 



3.7' 

(e,oo) 

K,oo) 



y = F 2 / 3 {x,y), y(0 = v 



and 77 



(0,0 
(o,£] 



Additionally, the graphs of solutions y Hmit contain the points (1,77£ 2 / 3 ). Hence there is no need to 
consider initial values different from y(l) = tj = yi Im i t (l) for 77 $J 0. 

Therefore all (derivatives y of) solutions to (fl) qualitatively have the shape of the dashed line. 
(y < 0, y' = /" ^ 0, furthermore / and /i imit have the same slope in 1). In (0,1), however, / 
decreases more rapidly than / limit , in (l,oo) less rapidly. 







The question now is: Are there other solutions of Note that e.g. the attempt to find a 
solution by solving y' = F(x, y) := F 2 /3(x, y)+s for some positive e leads to solutions that satisfy 
/ < in a bounded neighbourhood of £ only and not on the whole of R+. However, F(x,y) := 
F( 2 /3)+ a (x, y) = - |- (| + a) = -&&y for a > gives y a (x) = ry-exp (- /* ^±2 dt) = r,-x~(i +a ) 
as solution to 

/ 3a + 2 

Va = o Va and Vat 1 ) = V < 

3x 

and hence we obtain the following family of solutions to ([T]) : 

Lemma 2.4. (Family of solutions) For arbitrary c ^ 0, d € R, a £ [0, oo), the family of 
functions that is defined on H + by 

(d + c-sT>- a , for at [0,1/3) 
fa(s) ■= \ d + c ■ Ins , for a = 1/3 

[d~c-s^ a , for a € (1/3, oo) 

has the property that f a o dot: P§ym(3) — ¥ R, C <— > / a (detC) is convex. 

Remark 2.5. Although this condition is not necessary for the convexity of / o det, at least quali- 
tatively all solutions to ^ have the shape indicated in the graph below. (As we discussed in this 
section.) 



fa(s) 




























































































> 
















3IF 1 


"/ 3 -3 


















Ins 

/6 + 6 



/ 1Joiit (s) = -3sV3 + 3 

Acknowledgement 

We thank John Ball for his interest in this result. 



(i 



References 

[1] J.M. Ball, Constitutive inequalities and existence theorems in nonlinear elastostatics., Herriot 
Watt Symposion: Nonlinear Analysis and Mechanics. (R.J. Knops, ed.), vol. 1, Pitman, 
London, 1977, pp. 187-238. 

[2] , Convexity conditions and existence theorems in nonlinear elasticity., Arch. Rat. 

Mech. Anal. 63 (1977), 337-403. 

[3] , Some open problems in elasticity., Geometry, mechanics, and dynamics. (P. Newton 

et al., ed.), Springer, New- York, 2002, pp. 3-59. 

[4] J.R. Magnus and H. Neudecker, Matrix differential calculus with applications in statistics and 
econometrics, Wiley Series in Probability and Statistics, John Wiley & Sons, 1999. 

[5] L. Mirsky, An introduction to linear algebra, Clarendon Press, 1972. 

[6] P. Neff, 5 lectures, Advanced school on "Poly, quasi and rank one convexity in applied me- 
chanics"., Organizer: J. Schroder and P. Neff, Centre International des Sciences Mecanique 
(CISM), Udine (Italy), 24.-28. September 2007. 

[7] T. Rockafellar, Convex Analysis., Princeton University Press, Princeton, 1970. 

[8] J. Schroder and P. Neff, Invariant formulation of hyperelastic transverse isotropy based on 
poly convex free energy functions., Int. J. Solids Struct. 40 (2003), no. 2, 401-445. 

[9] J. Schroder, P. Neff, and V. Ebbing, Anisotropic polyconvex energies on the basis of crystal- 
lographic motivated structural tensors., J. Mech. Phys. Solids 56 (2008), no. 12, 3486-3506. 

[10] G. Strang, Inverse problems and derivatives of determinants, Archive for Rational Mechanics 
and Analysis 114 (1991), 255-265, 10.1007/BF00385971. 



7 



